Resonance scattering and singularities of the scattering function 
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Recent studies of transport phenomena with complex potentials are explained by generic square 
root singularities of spectrum and eigenfunctions of non-Hermitian Hamiltonians. Using a two 
channel problem we demonstrate that such singularities produce a significant effect upon the pole 
behaviour of the scattering matrix, and more significantly upon the associated residues. This mecha- 
nism explains why by proper choice of the system parameters the resonance cross section is increased 
drastically in one channel and suppressed in the other channel. 
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I. INTRODUCTION 

It is well established knowledge that the shape of mea- 
sured cross sections, in particular that of resonances, is 
described satisfactorily by poles of the scattering func- 
tion in the complex energy plane [l|. These poles corre- 
spond to particular solutions of the Schrodinger equation 
at complex energies by imposing the boundary conditions 
that the wave function be regular at the origin and have 
only outgoing waves at large distances @. They have 
been discussed recently again in the context of complex 
potentials Q. In fact, while for a hermitian Hamiltonian 
such solutions can occur only at complex energies with 
non- vanishing negative imaginary part [ij , a complex po- 
tential (associated with absorption) can give rise to such 
states at energies with vanishing or even positive imagi- 
nary part. As a rule, these specific solutions give rise to a 
pole of first order in the Green's function (E — ff) _1 and 
thus in the scattering function. This remains true also 
in multi-channel problems, even if there is a degeneracy, 
where the residue of the pole is a matrix of rank larger 
than one. 

There are, however, special singularities of the Green's 
function and thus of the scattering function being of a 
very different nature: the exceptional points (EP) of 
a Hamiltonian. They are also specific solutions of the 
Schrodinger equation or, in more general terms, of an 
eigenvalue problem, but for a non-hermitian problem. 
They have been discussed at great length theoretically 
0-01 and in a large variety of applications such as in 
atomic physics [1, |||, in optics [l(|, in nuclear physics 
and in different theoretical context in "PT-symmetric 
models [Hj], to name just a few. Depending on the partic- 
ular situation EPs can signal a phase transition [l3], [l4[ ■ 
The experimental verification of the existence of EPs re- 
quires a careful tuning of parameters of the particular 
open system investigated (see below Sect lll Bj ). Recently 



the topological structure of the square root branch point 
associated with an EP has been shown experimentally to 
be a physical reality [HtSt]]- 

While some of the quoted papers report about partic- 
ular physical effects of the EPs, either expected or mea- 
sured, no emphasis is placed upon the important role of 
the eigenstates when the vicinity of an EP is explored 
in the laboratory. We recall that an EP is characterised 
by the coalescence of not just two eigenvalues but also of 
their associated eigenstates; moreover, the norm of the 
associated eigenstate - there is only one - vanishes at the 
EP. It is the main purpose of the present paper to il- 
lustrate the dramatic effect of this vanishing norm upon 
the residues of the scattering function (which has a pole 
of second order at the EP [18|]) resulting in a distinctly 
different resonance behaviour in the different channels. 
To facilitate the demonstration we restrict ourselves to a 
two-channel problem. 



II. SHARP RESONANCES AND 
EXCEPTIONAL POINTS 

A. The Model 

To illuminate the basic mechanism of the sudden in- 
crease of a resonance cross section, associated with EPs, 
we strip the great variety of models down to the essentials 
and elaborate the universal cause by studying a generic 
model. It is well established that, in the close vicinity of 
an EP, a two-dimensional matrix model suffices to cap- 
ture all essential features associated with the singularity. 
We thus begin with the model Hamiltonian 



H(\) = H + H 1 (X) = H Q + XV 
wi \ , , / ei S 
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FIG. 1: Crossing of resonances without coupling (5 — 0). 
The three figures are taken for A = 0.53, KA„. oss = 0.5625 
and 0.59, respectively. The units are arbitrary but the same 
for all three figures. 
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where the parameters LUk and determine the non-inter- 
acting resonance energies = lj^ + Ae^, k = 1,2. They 
are chosen complex and such that the E/~ have negative 
imaginary parts; we focus on A between 0.5 and 0.6. As 
we arc interested in the effects of crossing and coalescence 
in scattering, we consider in the following the scattering 
matrix, viz. 

T i>k = #i(A) iifc + (Hi (A) (E - H(X))- 1 fli(A)) ilifc . (2) 

Since we mimic a multi-channel problem by an effective 
two-channel Hamiltonian, i.e. Ui and are complex, the 
scattering amplitude © does not satisfy the unitarity 
condition [l[. 

To describe, for (5 = 0, the simultaneous observation 
of the two resonances in both channels on equal footing, 
we have to consider the scattering matrix in a basis ro- 
tated by the angle ±7r/4 rather than in the basis given 
by (p}. That is, the labels i,k in the scattering matrix 
© refer to the basis {cl/V2 ± c2/V2, cl/%/2 =p c2/\/2} 
with {cl,c2} being an eigenvector of ([T]). This specific 
observational basis is essential to ensure that even for the 
non-interacting case (8 = 0) both resonances are equally 
and simultaneously present in either channel; the differ- 
ence between the two channels shows when 8 is switched 
on as is disccussed in the following sections. We choose 
the parameters and for convenient demonstration, 
that is wi = 1.55 - 0.007i, w 2 = 1.1 - 0.007i, ei = 
-0.4 - 0.0006i, e 2 = 0.4 + 0.0005i. The values have no 
particular significance, they serve to illustrate the prin- 
ciple, that is the effect of a near EP upon scattering; 
any other set that invokes crossing and coalescence of 
resonances would serve the same purpose as long as the 
imaginary parts of the poles are near to the real axis. 

B. The interacting resonances 

From © we obtain the cross sections ~ (TT^)kk as 
a function of the energy for various values of A. This is 
illustrated in Fig.l, where (TT^)n is plotted for 5 = 0. 
Without coupling, when the energies cross (a usual de- 
generacy), the figures for (TTt) 22 coincide with those 
of (TT Such crossing of resonances is often asso- 
ciated with a particular symmetry, this aspect is of no 
significance for our purpose. The crossing happens at 
A cross = 0.5625 - 0.00077i. 

The difference between (TTt)n and (TTt) 22 changes 
dramatically when the coupling 8 is turned on. In our 



FIG. 2: Energy eigenvalue trajectories when A sweeps from 
0.53 to 0.59 for 8 = 0.0115i. If 8 = the straight lines con- 
tinue as straight lines avoiding the peaks. Note that for our 
choice of parameters the trajectories run in opposite direc- 
tions, the respective starting points are indicated by a dot. 



example we have chosen a purely absorptive coupling 
as then the two EPs avoid the real A-axis (recall that 
the Hamiltonian is non-Hermitian even for 8 = ow- 
ing to the unperturbed widths; it means EPs may oc- 
cur for real A). For 8 = 0.0115i the EPs are at A + = 
0.579 - 0.00082i and A_ = 0.547 - 0.00073i with the 
energies 1.325 - 0.007029i and 1.325 - 0.007027i, respec- 
tively. The EPs sprout from A cross , a critical value used 
below is A cr = 5R(A + |A_)/2« A cross . The dependence 
on any parameter is highly sensitive when 8^0 owing 
to the proximity of the EPs invoking large derivatives of 
spectrum and eigenfunctions. In particular, the depen- 
dence on 8 is crucial and can change patterns as is further 
discussed below. 

In Fig. 2 we illustrate the trajectories of eigenvalues of 
the Hamiltonian (Q} obtained by variation of A for the 
parameter set considered above. We note the sharp de- 
viations from the straight lines of the trajectories bring- 
ing one resonance close to the real energy axis while the 
other is moving into the opposite direction. The two 
EPs lie basically within the centre of the two peaks; at 
each of these points the levels behave to lowest order as 
E± ~ Eep ± const. -\/A — Xep implying a large deriva- 
tive. Moreover, the state vectors have a similar strong 
dependence: while they are properly normalised their di- 
rections change rather swiftly when A is sweeping over the 
critical value A cr = 0.563. We emphasise that the phe- 
nomenon of width repulsion associated with level cross- 
ings as illustrated in Fig. 2 is usually not encountered for 
real values of the strength parameter A. It occurs in our 
model as both EPs are on the same side of the real A- 
axis being a consequence of dealing with a non-Hermitian 
Hamiltonian from the outset 1!) . The combined effect is 
the dramatic change for both cross sections in the imme- 
diate vicinity of the critical point as illustrated in Fig. 3. 

The obvious effect is the huge increase of the cross sec- 
tion for (TT^)n which is due to the narrow resonance. In 
contrast, there is a substantial drop for (TT')22 in com- 
parison with the centre figure of Fig.l. Moreover, the 
shape of this latter 'resonance' is no longer of the usual 
Lorentzian type (as was first noticed in Ref.jjoj). This 
dramatic change essentially takes place, albeit continu- 
ously, within the narrow interval {A cr — 0.03, A cr + 0.03}. 
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FIG. 3: Energy dependence of (TT f )n (left) and (TT t ) 22 
(right) for 8 = 0.0115i at the critical point X cr = 0.563. Note 
the huge difference in scale between the two drawings. The 
units are as in Fig.l. 



In fact, for the same values of the coupling the curves at 
A = 0.53 or A = 0.59 look just like those in Fig.l, where 
8 = is considered. This is understandable as the trajec- 
tories shown in Fig. 2 are virtually unaffected by the small 
coupling outside the interval {A cr — 0.03, A cr + 0.03}. The 
same holds for the eigenvectors. We note that the "peak 
on the peak" on the right hand side of Fig. 3 is again an 
effect of fine tuning. It is brought about by the conspir- 
acy of the two poles and their respective residues as they 
occur on the top and bottom cusp in Fig. 2. Decreas- 
ing for instance 8, the "peak on the peak" can disappear 
but the resulting shape still is not a Lorentzian, it rather 
retains the fiat plateau without the "second peak" . 

Of course, there is nothing special about the 11- 
channel that becomes so strongly enhanced. When the 
sign of the rotational angle 7r/4 is changed, the two eigen- 
functions simply swap and it is then the 22-channel that 
appears so dramatically enhanced while the 11-channel is 
suppressed; similarly, as the phases of the (normalised) 
eigenvectors of ([T]) are arbitrary, changing the sign of 
either cl or c2 will also swap the roles of the two chan- 
nels. The important point is the strong enhancement 
of one channel while the other is suppressed. In either 
case, the effect occurs only in the immediate vicinity of 
A cr = 0.563. The pattern in Fig. 3 has its origin not only 
in the square root branch points of the spectrum but 
rather in the associated eigenstates. This aspect, that 
is the effect of the eigenstates upon the residues of the 
Green's function, is often overlooked in the literature; it 
is not only the pole positions but also the strong and dis- 
tinctly different variation of the two residues that brings 
about the effects discussed here. 

The trajectories of the spectrum in the energy plane as 
illustrated in Fig. 2 provide the major clue for an under- 
standing. Within the range of interest for real A-values 
the two peaks pointing in opposite directions are well 
described by 



E±(\) nE r ± V(A-A_)(A- A+), 



(3) 



where 



A± = T: 



llUJi - LU 2 ) 



2S±i(e 1 -e 2 ) 
are the two EPs. For small 8 they lie near to each other 



and degenerate into a single diabolic point when the cou- 
pling parameter 8 vanishes. The constant E r ensures that 
both energies have a negative imaginary part for real A 
to allow their interpretation as true resonances. Notice 
that care has to be taken for the values to remain in the 
same Riemann sheet in a numerical generation of a plot 
like the one in Fig. 2 where real values of A are used. At 
the critical point X cr - that is where the two tips of the 
two peaks occur - one energy value is very close to the 
real axis while the other is more remote from the real 
energy axis. These two pole terms occur in the Green's 
function, and it is here where the corresponding residues 
play a decisive role. The analyti cal expressions for the 
residues blow at the EPs like 1 / \/(A — A+)(A — A_) ow- 
ing to the vanishing norms of the eigenstates. We give a 
simplified expression representing numerically the correct 
analytical terms quite well for real A. It reads 



r±( A)«l±— A 

2 ^A-(A+ + A_)/2 



(4) 



where the constant A is suitably adjusted yielding 
|r+(A)| < 1 at A cr . A combined diagram of the real func- 
tions |r+(A)| and |r_(A)|, as a function of 5iA, looks sim- 
ilar in shape to Fig. 2 with the difference that the peaks 
reach near unity and zero, respectively. The diagonal 
elements of the full Green's function are then given by 



{E-H)ll = 
(E-H)£ = 



r+(A) 



r_(A) 



E-E + (X) E-E-(X) 

r-W r + (A) 
E-E + (X) E-E-(X) 



(5) 



while the correspo nding poles of the off-diagonal terms 
have the residues -y/r+(A) • rZ(Xj. With the insight from 
Eqs. (|3I4I5|) we qualitatively understand the shapes of 
Fig. 3. The strong peak (left) originates from the near 
pole at .E+(A cr ) with its strong residue r+(A). The right 
hand figure is produced by the remote pole at B_(A r ) 
with its strong residue while the small peak on top is 
generated by the near pole with the very small residue. 
The rather dramatic A— dependence is due to the prox- 
imity of the EPs; of course this is related to the actual 
value of the coupling 8 and will have a sensitive depen- 
dence on 8 at A = A cr . Note that, if the coupling in 
our two-dimensional model is further increased, the nar- 
row resonance is eventually crossing the real energy axis 
and acquires a positive imaginary part. It then can no 
longer be interpreted as a decaying state, that is as a 
physical resonance. It may be interpreted, for example, 
as an inversely populated level producing light amplifi- 
cation generated by a radiation process in laser physics 

m. 



III. CONCLUSION 

Let us discuss a few cases where the phenomena de- 
scribed in this paper have been encountered. In a model 
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for quantum transport p2j |. transmission increases dras- 
tically for a specific width of some external barriers simu- 
lating the coupling strength of the system with the envi- 
ronment; this parameter plays the role of our A. The au- 
thors stress that the complex eigenvalues of the effective 
Hamiltonian are very important for understanding the 
transport properties of the system. The analysis is based 
in part on findings of [23| where the sensitive dependence 
upon the threshold energy of the continuum coupling for 
loosely bound nuclei has been demonstrated using an ef- 
fective two by two non-Hermitian Hamiltonian. In Q a 
specific frequency of the laser beam leads to the occur- 
rence of a 'resonance' pole with zero width; within its 
immediate vicinity this parameter can be rewritten as 
our parameter A. The quoted paper concentrates on spe- 
cific solutions of the Schrodinger equation, which are in 
fact the Gamow states (only outgoing waves as boundary 
conditions). They correspond to the poles of the Green's 
function in our two-level model (which has of course no 
continuous spectrum). In all these cases, the extremely 
strong dependence on the coupling of the pole positions 
close to or upon the real axis, is the effect of near excep- 



tional points. The openness of the systems is crucial for 
all these phenomena. 

In summary, the dramatic features of particular reso- 
nance behaviour are explained by the square root singu- 
larities of spectrum and eigenstates, especially the van- 
ishing norm of the latter that gives rise to the dramatic 
and distinctly different behaviour in each channel. In 
a two channel problem, that was chosen as a typical 
case, they invoke strong dependence on the interaction 
strength and significant deviations from the usual pat- 
terns associated with isolated or overlapping resonances. 
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